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1. Introduction 

Lifshitz points (LP) are familiar examples of multi-critical points. At a LP a disordered, 
a homogenous ordered, and a modulated ordered phases meet [1, 2, 3, 4]. In the case of 
systems with m-axial LP, there is a degeneracy such that long-range order modulated 
along any of m distinct axes can occur in the modulated ordered phases. The study of 
critical behavior at such LP began immediately after their discovery in the middle of 
the 1970s. 

Unfortunately, the technical difficulties one is faced with in analytical 
renormalization group (RG) calculations are enormous. This is the main reason why RG 
results based on systematic expansions, such as expansions in powers of e = d*{m) — d 
about the upper critical dimension d*{m) = 4 + m,/2, or in powers of 1/n, where n 
is the number of components of the order parameter, had remained quite scarce for 
decades. Furthermore, early e-expansion results obtained by two different groups (cf. 
[5] and [6, 7]) had yielded contradictory results, and these discrepancies had remained 
unclarified for many years. The results of reference [5], which were restricted to the 
special cases of bi- and hexa-axial LP m = 2 and m = 6, were reproduced twenty years 
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later by field-theoretic means [8]. However, a full two- loop RG analysis in d*{m) — e 
dimensions was reported only in 2001 [9, 10]. This gave the e expansions to second 
order of all four main independent critical exponents as well as the correction-to-scaling 
exponent for general values of m, besides resolving the mentioned discrepancies. || 

In a recent paper [13] (hereafter referred to as I), we have shown how the 1/n 
expansion can be applied to the study of critical behavior at m-axial LP. We determined 
the correlation exponents 77^,2 and tjla, and the related anisotropy 9 of d-dimensional 
systems to first order in 1/n for general values of m. The results took the form of 
complicated integrals whose integrands involve further multi-dimensional integrals. We 
were able to check that in the isotropic limits m ^ and m ^ d they correctly reduce 
to known results, namely, the expansions to order 1/n of the Fisher exponent r/ at a 
usual critical point [14, 15], and that of 77^4 at the isotropic rf-axial LP [16], respectively. 
Furthermore, we could show that our 0{l/n) results are in conformity with published 
dimensionality expansion results about the lower critical dimension d^:{m) = 2 + m/2 
both for m = [17] and for m = 2 and m — 6 [18]. 

The virtues of the 1/n expansion are well known: It can be applied in arbitrary 
fixed dimensions d, does not rely on the smallness of a further expansion parameter 
such as e, and yields nontrivial results below the upper critical dimension d*{m) 
in a mathematically controlled fashion. Besides its capability of providing valuable 
information about the critical behavior for given d and m, it allows for nontrivial checks 
on e-expansion results. Unfortunately, the complicated general form of our results in I 
prevented us from proving their consistency with the e-expansion results of [9, 10] for 
general m. We could verify it in the isotropic limits m — > and m ^ d. However, the 
only other case in which we could explicitly demonstrate this consistency by analytical 
means was that of m = 2. This is unfortunate for at least two reasons: first, it excludes, 
in particular, the physically important uniaxial case m — 1; second, unlike both the 
contributions of first order in e = d*{m) — d as well as their 0{d — d^{m.)) analogues 
in the expansions about the lower critical dimensions d^:{m), the O(e^) terms exhibit 
a nontrivial m-dependence [9, 10]. This propagates into the 0(l/n) series coefficients, 
and ought to be checked. 

The purpose of this article is to fill this gap and prove that our large-n expansion 
results [13] to order 1/n for the exponents riL2{d,m,n), riLA{d,m,n), and 9{d,m,n) are 
fully consistent with the e-expansion results of references [9, 10]. In the next section 
we first provide the necessary background, recalling the continuum model on which our 
analysis is based as well as our results for general m given in I. We then show that these 
results can be rewritten in a form allowing analytic comparisons with the O(e^) results 
of references [9, 10] for general m. The proof that they are in conformity with the latter 
is given in section 3. The closing section 4 contains a brief discussion and concluding 
remarks. 

II Alternative results reported in reference [11] could be refuted [12]. 
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2. Lcirge-n expansions of the correlation exponents at the Lifshitz point 

Just as in I, we consider a model defined by the Euclidean action 



7i [</,] = 1 j d'^-"^r J (Tz 



(1) 



Since we intend to work directly at the LP, we have set the coefficients of the quadratic 
terms to their corresponding critical values tlp and plp If- 

Here </> = (t>{x) is the usual n-component order-parameter field. Its d-dimensional 
position vector x — (r, z) e has a (d — m) -dimensional "perpendicular" component 
r and an m-dimensional "parallel" one, z. The subspace associated with z is the 
one in which modulated order can occur in the corresponding phase; that of r is 
its orthogonal complement. A similar decomposition has been made for the gradient 
operator V = (Vr, V^). Thus is the Laplacian in the parallel subspace. 

Employing the notational conventions of I, we write the wave- vector conjugate to 
X = {r,z) as k = {p,q), with p G M'^^™ and q G M"'. Further, we introduce the 
two-point cumulant G^{r, z) and its Fourier transform G^{j), q) through 

G^(r, z)^- mr, z) ■ 0(0, 0)) - z)) ■ (0(0, 0))] 

77" 

= / / e^(^-P+^-^)G^(p,g), (2) 

J p J q 

where J^'^"™^ = (27r) /^^.^ and j^^'^ = {2%)-"' J^^d'^q denote normalized 

{d — m)- and m-dimensional integrals, respectively. 

As discussed in I, the full propagator (5^(p, q) becomes a generalized homogeneous 
function G^^\p, q) in the hmit of large length-scales. The latter function satisfies at the 
LP the homogeneity relations 

gJ^Hp, q) = ^(1, qp-') = ^(-^(pg-V^, i) . (3) 

Only two exponents are independent here since the usual scaling relation 

e = l^ (4) 

must hold for the anisotropy index 9 by consistency. 

In the limit 7i — > 00 with nX = fixed, G^jf^ij), q) reduces to the Gaussian propagator 
&'^\p, q) pertaining to the Hamiltonian (1) with A = tlp = Plp = 0. We have 

lim Gf\p,q)^&^\p,q)^-^,. (5) 

nA=const " " 

% In the following it is tacitly understood that d with d < d*{m) = 4 + m/2 and m are chosen such 
that a LP exists. This requires, in particular, that d exceeds the dimension 2 + to/2 below which the 
homogeneous ordered phase becomes thermally unstable because of spin-wave excitations. However, it 
also requires that the modulated ordered phase remains stable. For a discussion of these delicate issues, 
see the review article [4] and its references. 
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At order 1/n, self-consistent equations must be solved which were discussed in I and 
need not be repeated here. To this end, we looked for solutions of the scaling form (3), 
utilizing the ansatzes 

VL2 = ^ + 0in-'), VL, = ^ + 0{n-'), e^U — + 0{n-'), (6) 
n n 2 n 

together with corresponding 1/n expansions for the scaling functions in equation (3) 

and the relation 

(1) (1) 

implied by the scaling law (3). This led to consistency conditions (equations (27) and 
(28) of I) from which we obtained the results 

(1) Ka-m A^™^ 2V,{q') 

r 

and 



(1) /-^'-^^ 8V2{p') 1 



"^^^ 4m(m + 2) {p^ + 1)5 I{p, 1) ' 

Here Kd-m and Km, defined by 

""^wro with (10) 

where So is the area of a unit sphere in D dimensions, are conventional factors resulting 
from the angular integrations ect D — d — m and m. Further, Vi{q^) and V2{p^) denote 
the polynomials 

Vi{q*)=4:-{d-m){l + q^) (11) 

and 

V2{p^) = 3(8 - m)(6 - m) + 5(m^ + 2m - 96)p^ 

+ (m^ + 50m + 144)p^ - m(m + 2)p^ . (12) 

Finally, /(p, q) represents the analogue of Ma's "elementary bubble" XIHX [14] for 
the LP: 

Md-m) Am) ^ 

whose homogeneity property 

I{p,q) =p-'I{l,qp-'^^) =q-^'I{pq-\l) , e = 4-d + m/2 , (14) 

we recall for later use. 

Let us first show that the results (8) and (9) for the 0(l/n) coefficients can be 
rewritten as 

""^^ - 2{d -m)L p^ + 1 I{p, 1) ^''^ 
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and 

(1) ^ Kd^rn f'-'^^ 1 4 1 , . 

"^^^ 2m(m + 2) 1 + " I{l,q) ' ^ ^ 

respectively, where Vg = (Vg)^. These representations of the coefficients r)^j^2 VlI 
are well suited for determining the 0(e^/n) contributions to the exponents rjL2 and 77^,4. 
They will be employed as starting point in our proof of consistency with the e-expansion 
results of [9, 10] given in the next section. 

To derive these representations, note first that the actions of the D-dimensional 
Laplacian = d^/dK^ and its square on functions f{K^) of i^^ = ^^^^ 

can be written as 

V]^f{K^) = 2Df\K^) + AK^f\K^) (17) 

and 

Vy{K^) = 4L>(L> + 2)f"{K^) + 16{D + 2)^^^^^ (^') + IQK^f^^ (K^) , (18) 

where f^^\.) means the sth derivative of the function /(.). 

Using these relations, it is straightforward to see that the rational functions 
appearing in the integrands of (8) and (9) can be expressed as 



(1 + g4)3 p2^q4 
and 

8^2(P^) ^4 1 



(19) 

P2=l 



V 



[1+p^y "^P^ + Q^ 



(20) 

We now insert these results into equations (8) and (9), use hyper-spherical 
coordinates for the integrals fj^"^^ and f^^ make the changes of variables q 
p = and p q = p^^^"^ in the radial integrals over q and p, and utilize the 
scaling property (14) to express /(l,p~^/^) and I{q~^, 1) in terms of I{p, 1) and /(I, q) 
respectively The derivative term on the right-hand side of equation (19) becomes 
p'^Vp {P'^p^ + l)"^|p=i = Vp {p^ + 1)~^. The Vq term in equation (20) transforms in 
a corresponding fashion. One thus arrives at expressions that agree with equations (15) 
and (16) except that the derivatives act to the left. Integration by parts then yields the 
claimed results. 

A straightforward, though important first application of them is to show that the 
e expansion of the coefficients 'q^H and -q^H starts at order e^: 

^S,4 = + 0(6^) . (21) 

To see this, note that I{p, q) has a Laurent expansion about e = of the form 

/(p,5) = ^ + /o(p,g) + 0(e) (22) 
with a momentum- independent residuum given by 
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according to references [9, 10] (see equations (7), (24) and (89) of [9] or (38) and (39) 
of [10], where /_i was denoted Fjn,o)- Hence 



72 1 _ 2'^P^0(p,g) 



with Vp/o(p, q) = [Vp/(p, ?)]e=o- Analogous results with replaced by hold. Thus 
both ?7^2^ and 77^4 are indeed of order e^, and for their O(e^) expansion coefficients 77^2,4 
we obtain from equations (15), (16) and (24) the results 

and 



'fif^-^^^^-k- r'r^vl/(i,,)U. (26) 



2m(m + 2) P^Jq 1 + ' ^ " 
3. Epsilon expansions of r]^^2 VlI 

We are now ready to present the announced proof of consistency We shall show that 
the e expansions of the correlation exponents 77^,2 and 77^,4 are related to the coefficients 
^2 and 77^4^ via 

^^2,4 = ^^^r;S,4 + 0(^')- (27) 

Since the prefactor (7i + 2)/(7i + 8)^ on the right-hand side reduces to I/ti in the large-7i 
limit, consistency between the e expansions to second order and the 1/n expansions to 
first order is an immediate consequence. 

The O(e^) results of references [9, 10] for the exponents ?7l2,4 involved single 
integrals, which for general m had to be computed by numerical means. In the notation 
of the second of these publications, they read 



and 



where 



and 



j^{m) = Bm dv v"^-' ^^{v; m, d*) , (28) 
Jo 

j,{m) = B^ dv ^%v; m, d*) , (29) 
Jo 



= ^±1^ (30) 
-'-1 



^v;m,d) = G^'^\l,v) . (31) 

The latter is the scahng function associated with the free position-space propagator 
G^^\r,z), whose scaling properties 
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we recall. For general values of m and d, it is a difference of two generalized 
hypergcomctric functions 1F2. This is why no analytic results for the integrals jtf,{m) 
and ja{ni) are available for general m. In reference [10], the reader may find numerical 
results for them at m = 1, 2, . . . , 7 along with analytical ones for m = 2 and m = 6. 
To prove the relations (27), we must show that 

V^l?-^3M) (33) 



m 



and 

^m'^ = 7^ 7 jM) ■ (34) 

'^^ 2m(m + 2)-^''^ ' ^ ' 

Let us start from equation (25). Its integral J^"^ has the form of a scalar product 
{j\g) in L2(]R^~™'^^), the space of square integrable functions, that is evaluated in the 
p-representation. In r-space the bra (/| is represented by /(r)* = where /(r) is 

the Fourier qr-transform of G'*^'^-'(r, z\ taken at an arbitrary unit qr- vector q. Performing 
the angular integrations in the required m-dimensional integral is straightforward and 
yields 

//•oo 
d^zG^''\r, z) e'^-^ = (27r)"^/2 j dz z"^/^ Jf-i{z) G^^'^r, z) ,(35) 

where from now on e is set to zero in G^^\r, z). 

Likewise, {r\g) = g{r) is the Fourier q-transform of the function (— r^) [G^^\r, z)]"^ 
for q — q. In the resulting expression for 

/POO POO 
^4-m/2^ dz dz'r^ZzT^^ Jin_^{z)G^''\7 

Jo Jo ^ 

xJ^_,{z')[G^'\r,z')r (36) 

substitute the first of the scaling forms (32) along with equation (31). We then make 
the changes of variables r ^ v = zr~^^'^ in the radial part of the integration over r and 
z' ^ ( = z'/z. The resulting integral over z is the special case of the closure relation 
for Bessel functions^ 

POO 

/ dzzJ,{Cz)J,{bz)^S{C-b)/C (37) 

with n = m/2 — 1 and b — 1. The integral over C, can now be performed. Substituting 
the result into equation (25) and noting (28) and (30) then gives the asserted result (33) 
for the coefficient rj^^f^ . 

The corresponding expression (34) for r])^^ can be proven in an analogous fashion. 
The integral j''^^ is a scalar product {h\w) in L2(M"') between h{z), the Fourier p- 
transform of G{r^z)^ and w{z)., that of [G'(°)(r,^)]2, taken at a unit p-vector p. We 
perform the angular integrals in the Fourier integrals J (i'^^™/^r and J (i'^~™/^r', and 
make the changes of variables z ^ v = zr~^l'^ and r' ^ ^ = r'/r. The integral over r is 
of the form (37) with = 1 — m/4 and 6=1. Once the integral over C is performed, 
the desired result follows from equations (26), (30), and (29). 

See, for instance, [19] 
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4. Concluding remarks 

In this paper we have shown that the large-n expansion yields 0(l/n) results for the 
correlation exponents r]L2 and r^L4 and the related anisotropy exponent 9 for general m, 
which are fully consistent with the e-expansion of references [9, 10]. In view of the long- 
standing discrepancies mentioned in the Introduction and the great technical challenges 
encountered in both expansion methods beyond lowest order, the established consistency 
is very gratifying, providing nontrivial checks of the results of both expansions given in 
references [9, 10] and [13], respectively. 
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